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Abstract 

The observed microwave background anisotropics in combination 
with the theory of quantum mechanicahy generated cosmological per- 
turbations predict a well measurable amount of relic gravitational 
waves in the frequency intervals tested by LISA and ground based 
laser interferometers. 



*Based on a talk presented at the First International LISA Symposium, 9-12 July 1996, 
Rutherford Appleton Laboratory, Oxfordshire, UK 
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At the first glance, the ground-based and space laser interferometers 
for gravity wave observations, as well as the Weber bar technique, do not 
have much in common with the ongoing and planned radio-astronomical 
measurements of the microwave background anisotropics . The solid-state 
detectors are sensitive to gravitational waves in the 10^ Hz frequency range, 
the laser interferometers are sensitive, correspondingly, to the frequencies 
(10 — 10^) Hz and (10~^ — 10"^) Hz, whereas the microwave background 
anisotropics directly reflect only the variations in the cosmic temperature 
and, if they are caused by gravitational waves, can only provide us with 
information about extremely low-frequency gravitational waves - (10^^^ — 
10^^*^) Hz and lower f^. However, the basic physics which enables us to 
see gravitational waves with the help of laser interferometers or microwave 
background anistropies is exactly the same: alterations in frequency and 
phase of an electromagnetic signal propagating in the field of gravitational 
waves. Most importantly, the relic stochastic background of gravitational 
waves, to be discussed below, extends from very high frequencies of the order 
of 10^ Hz to extremely low frequencies of the order of 10~^^ Hz and lower. 
It can be observed by all these techniques, and the predictions about the 
expected gravitational wave amplitudes in various frequency intervals are 
connected to each other by the theory [^]. 

Although the measured large-angular-scale anisotropics in the cosmic mi- 
crowave background radiation (CMBR) [0] could be expected on general 
grounds, their actual existence raises certain theoretical problems. The ob- 
served Universe is far from being homogeneous and isotropic, but becomes 
more and more so when one expands the study to larger and larger scales. 
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The anisotropies signify the presence in the Universe of cosmological per- 
turbations with very small amplitudes (the dimensionless deviations are of 
the order of 10~^) but with extremely long wavelengths, of the order of and 
longer than the present-day Hubble radius In, Ir ^ '^^ 10^*^ cm. Regardless 
of nature of the perturbations responsible for the observed large- angular- 
scale anisotropy, that is, regardless of whether they are mostly density per- 
turbations, or rotational perturbations, or gravitational waves, or all mixed 
together, there exists a puzzling question of their origin. The first wonder is 
whether they are remnants of the originally inhomogeneous and anisotropic 
Universe or, alternatively, were generated by some mechanism in the orig- 
inally homogeneous and isotropic Universe. Since the perturbations of our 
interest are weak, we can use the linearized Einstein equations for the de- 
scription of their evolution. 

It is difficult to maintain that these perturbations are simply survived. 
The photons of the CMBR have become free and started their journey to us 
sometime at the beginning of the matter-dominated era. In the preceeding 
radiation-dominated era, the general solution for a Fourier component of the 
metric perturbations is 



where A and x ^'^^ arbitrary constants, a(?7) is the scale factor of a FLRW 
universe 



and a{ri) ^ t] in the radiation-dominated era. Let us take the amplitude A at 
the level A ^ 10"^, in rough agreement with observations, and return back in 
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depth of the radiation- dominated era by sending rj to zero. The scale factor 
0(77) diminishes by at least the factor 10^ by the time of reaching the era of 
the primordial nucleosynthesis. The dangerous term A cos nrj sin xe'" ''/ 0(77) 
must be cancelled out with a great precision, if we do not want the hn{r], x) 
to become of the order of 1 and destroy in this way the homogeneity and 
isotropy of that era. That is, the phase x must be or tt with a precision of 
10~^, or even much higher if we want to proceed further back in time. The 
neighboring solutions must all have the same (or differing by tt) phases, that 
is, the distribution of the phases must be very narrow (highly "squeezed"); 
and the waves must be standing, and not traveling. (We will see below 
that all that is automatically guaranteed if the perturbations are generated 
quantum mechanically.) The solution (1) is strictly valid for gravitational 
waves, but the same argument is applicable for density perurbations as well. 

It is still possible that the perturbations of our interest are classical, 
deterministic remnants of a strongly inhomogeneous anisotropic universe of 
a very distant past. It is also possible that, for some miraculous reason, 
the phases have been chosen rightly with enormous precision, so that the 
perturbations are classical, deterministic remnants of a universe which was 
almost homogeneous and isotropic from the very beginning. To the present 
author, these possibilities do not seem to be likely, even if they can be shown 
to be consistent with all available data. We need to turn to posibilities of 
generating the perturbations in an originally FLRW universe. 

In principle, there are several options to do that. For instance, one can 
try to exploit the fact that the number of unknown functions of time partici- 
pating in the perturbed Einstein equations is always greater than the number 
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of equations. These functions describe nonadiabatic pressure, entropy per- 
turbations, anisotropic stresses, etc.. By manipulating with these functions 
and making additional assumptions, such as that these functions represent 
"cosmological defects" or "causal seeds", one can produce the required per- 
turbations, but essentially "by hand" . It seems to the writer that we should 
first try to build the theory on a minimal number of hypotheses. 

It is known that the Einstein equations plus basic principles of quantum 
field theory allow (in fact, demand) the quantum- mechanical generation of 
cosmological perturbations from the vacuum state, as a result of parametric 
(superadiabatic) interaction of the quantized perturbations with strong vari- 
able gravitational field of the very early Universe. This process is possible 
for gravitational waves 0, for density perturbations |^, §], and for rotational 
perturbations . However, in the last two cases we need to assume that the 
primeval matter was capable of supportng the oscillations and that they were 
properly coupled (similar to gravitational waves) to the "pumping" gravita- 
tional field. If, as is assumed in the inflationary hypothesis, the very early 
Universe was governed by a scalar field, and if the field was minimally cou- 
pled to gravity, the quantum mechanical generation of density perturbations, 
in addition to the inevitable generation of gravitational waves, was possible. 
Below, we will follow the line of the quantum-mechanical generation of cos- 
mological perturbations. 

In the presence of perturbations, the metric tensor and the energy-momentum 
tensor can be written in the form 
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(3) 




In cosmology, one usually considers the for simple models of matter, 
such as perfect fluids or scalar fields. The perturbations hij and Tj^^^^ are 
linked together by a set of the linearized Einstein equations. It is convenient 
to expand the perturbations over spatial harmonics e*""^, e~*""^, where 
n = {n^, n^,n^) is arbitrary wavevector and the wavenumber n is 



Demonstrating certain mathematical skill and understanding of the phys- 
ical side of the problem, one can show that, for the simple models of matter, 
the perturbed Einstein equations for each n-mode and for each of the three 
types of cosmological perturbations (density perturbations, rotational pertur- 
bations, gravitational waves) can be reduced to a single differential equation 
of second order. It is only this one equation that defines the dynamical con- 
tent of the problem and needs to be solved. All the components /iy , T^^^^ can 
then be found from its solutions by algebraic and differentiation/integration 
operations. This equation has the form of the equation for an oscillator with 
a variable frequency. The frequency varies due to the presence of the time- 
dependent scale factor a{r]) which plays the role of the gravitational pump 
field. The very form of this equation explains why the cosmological pertur- 
bations can be parametrically amplified (if the initial classical amplitude was 
not zero) or quantum-mechanically generated from the zero- point quantum 
oscillations. In case of gravitational waves, this equation is 
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fil + [n^ ]fin = (4) 

Oj 

and hn{r],x) ~ ^/inC*"-''. 

We will briefly summarize the main points of the quantum- mechanical 
generation of cosmological perturbations (for more details, see a recent paper 
|T0[| and references therein). The quantum- mechanical operator for hij{ri,x.) 
can be written as 



1 /"OO ^1 t 

^'^^'^'''^ = ^)J2^^ |^rf3n5:P,,(n)-^[c„(r^)e— + c:(r^)e— ]. (5) 

Each of the three types of perturbations has two polarisation states 
{s = 1,2) described by two polarisation tensors Pij. The normalization con- 
stant C is, up to a numerical factor slightly different for each type of the 
perturbations, the Planck length Ipi, Ipi = {Gh/c^Y^'^. For gravitational 
waves, C = y/lGirlpi. The time dependent creation and annihilation opera- 
tors c^irj), Cnirf) are governed by the Heisenberg equations of motion. The 
parametric nature of the interaction Hamiltonian allows one to apply the 
Bogoliubov transformation and to express the c^{'r]), Cn{ri) in terms of their 
initial values c^{0), Cn(0). The operators c„(0), Cn(0) define the initial vac- 
uum state |0) for each n and s: Cn|0) = 0. The field operator (5) acquires 
the form 



/^,,(r/,x) = |'^rf3^^^,^.(n)-^[L(r/)e— c„(0) + £(r?)e— cl(0)] 

(6) 



7 



where hn{v) essentially the solutions to the single second- order differen- 
tial equation (for each type of perturbations) mentioned above. 

In the Schrodinger picture, the initial vacuum state |On)|0_n), for every 
pair n, — n of modes, evolves into a two-mode squeezed vacuum quantum 
state. The modes affected by the amplification process will have a large 
mean value of the number operator N (large occupation number) and a large 
variance of N. The conjugate variance of phase will be highly squeezed near 
the values x = 0, tt (in the representation (1)). Classically, the generated 
perturbations can be treated as a stochastic collection of standing waves. 

The scale factor a{ri) of the very early Universe (well before the era of 
primordial nucleosynthesis) is not known. It depends on the unknown equa- 
tion of state of the extremely dense matter (we are not quite sure even about 
the equation of state in cores of neutron stars). However, it is likely that the 
evolution was significantly different from the law of expansion of a radiation- 
dominated universe. If so, gravitational waves must have been generated, and 
density perturbations, as well as rotational perturbations, could be generated 
too. 

Certain properties of the generated perturbations are universal, indepen- 
dent of a concrete form of a{rj) in the very early Universe, which we still keep 
quite arbitrary. These properties are determined by the fact that the per- 
turbations are placed in the squeezed vacuum quantum states. For instance, 
the expected (mean) quantum mechanical value of hij{'r),:x.) is zero in every 
spatial point and at every moment of time: {0\hij{r], x)|0) = 0. However, the 
variance is not zero and does depend on time: 
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(0|%(ry,x)/i^^(,7,x)|0) = nY^\hn{ri)?dn. (7) 

Eq. (7) defines the power spectrum 

nn)^^,nj2\Uri)?- (8) 

For numerical estimates, it is convenient to use the characteristic amphtude 
h{n) defining this amphtude as the standard deviation (square root of vari- 
ance) per logarithmic frequency interval: 



h{n) 



^2 2 ^ nV2 



s=l 



(9) 



It follows from Eq. (8) that for perturbations with wavelengths shorter 
than the Hubble radius, the power spectrum P{n) is not a smooth but an 
oscillating function of the frequency (wave number) n. Specifically for short 
gravitational waves, the generated stochastic background of the waves is 
not a stationary but a nonstationary noise, in the sense that the temporal 
correlation function of the field should depend on individual moments of 
time, and not only on the time difference. The large variance of N and hence 
the large variance of the amplitude of the perturbations will lead to large 
variations in the angular correlation function for the microwave background 
anisotropies, and so on. However, in this presentation, we will not go into 
the details of statistics and modulated spectra. We will be happy to show 
that we can get a right numerical level of the expected signal. 

The amplitudes and spectra of the generated perturbations depend on the 
strength and variability of the pump field, which in our case is completely 

9 



determined by the scale factor 0(77). Having the freedom of playing with the 
unknown part of 0(77) describing the very early Universe, one can derive, with 
the degree of completeness allowed by quantum theory, all the characteris- 
tics of the perturbations in the present Universe. For instance, some of the 
gravity-wave spectra derived in this manner are shown in Fig. 1 (adopted 
from m and updated). The original graph used, in the radio-astronomical 
fashion, the spectral flux density (ergs / sec cm^ Hz ster) as its vertical 
axis. In Fig. 1, we have used the vertical axis in terms of the dimensionless 
characteristic amplitude h{v), Eq. (9), where the dimensionless frequency n 
has been translated into the present day frequency v measured in Hz. The 
spectral cosmological fi-parameter due to the contribution of gravitational 
waves, fig(z^), is defined as [|] 




where is the critical cosmological energy density and vh ~ 10~^^ Hz is the 
Hubble frequency. Because the graph is the same, the vertical axis is now 
not universally homogeneous. 

In Fig. 1 one can also see the observational upper limits for stochastic 
gravitational waves, marked by arrows, that were valid at that time, 10 
years ago, as well as the expected sensitivities of then proposed new gravity- 
wave detectors. Some of the sensitivity curves, in particular for the space 
interferometer now called LISA, should be significantly modified, but we leave 
them in the old shape because the modifications will not be very important 
for our further discussion. 
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Figure 1 ; Possibk spectra of reik gravstaliona] waves. Thti soiid Ittie i!i [jie pTfttittitin 
based on the measured microwave, backgniuntf anisofjiofiies 



Let us start from the two spectra with the maxima near 10~* Hz and 
10^ Hz. These spectra were derived from bizarre cosmological models 
designed specially to produce the maximum possible amount of gravitational 
waves in the frequency intervals where the millisecond pulsar (MSP) and 
bar- detector techniques, then most favorite, were operating. One can also 
make equally bizarre assumptions about localized sources (such as colliding 
bubbles, phase transitions, string loops, strings with attached monopoles at 
the ends, etc.) and produce maxima virtually in any frequency interval, for 
the benefit of every individual experimental group. Generally speaking, we 
should keep eyes open to all these possibilities, they all are not forbidden. 
A different question is whether we will be surprised if the predicted signal is 
not detected and what we will learn from that fact. 

As was explained above, the quantum-mechanical (parametric) gener- 
ating mechanism relies only on the validity of general relativity and basic 
principles of quantum field theory. The law of expansion of the very early 
Universe is not known but this is what we will learn, or at least will place re- 
strictions on, by detecting or not detecting the predicted signal. For instance, 
the once popular cosmological model governed by matter with the stiff equa- 
tion of state p = e can already be ruled out, because the amount of the 
created high-frequency gravitons would be too big and would be inconsistent 
with available cosmological data ^ pT] . 

A large class of expanding cosmological models is described by the scale 
factor 

a(r7)=/o|r?r+^ (11) 
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where Iq and (3 are arbitrary constants. The rj time grows from — oo and (3 < 
— 1 at the initial stage of expansion. The constant /q has the dimensionahty of 
length and is effectively responsible for the Hubble-radius (time-dependent, 
unless (3 = —2) of the very early Universe. With this scale factor, the Einstein 
equations require the equation of state to be in the form 

For f3 = —2 one has p = — e, see Eq. (12), the case called inflation. 

Solving the gravity-wave equation (4) one can show that today's values 
of the characteristic amplitude /i(z/) should be as follows (ignoring the mod- 
ulation of the spectrum which takes place for z/ ^ u^)- 

For u < uh, 



- ■ (13) 

For i^// < z/ < z/m, where z/^ is determined by the time of transition from 
the radiation-dominated era to the matter- dominated era, ~ 10^^^ Hz, 



For t'm < < '^c, where Uc labels the highest frequency waves marginally 
affected by the amplification process and above which the spectrum sharply 
falls down, Uc ~ 10^ Hz in currently discussed models. 



M^)-t^ - - (15) 

The scale factors (11) which are power-law dependent on rj time generate 
spectra which are power-law dependent on frequency z/ p. 
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As we see, the numerical level of the predicted amplitudes depends on 
the fundamental constants G, c, h combined in Ipi and a couple of unknown 
cosmological parameters, such as /q and p. Of course, it is not for the first 
time that something observable is beautifully expressed in terms of the funda- 
mental constants and a couple of parameters only. For instance, the maximal 
masses and radii of white dwarfs and neutron stars depend essentially only 
on G, c, h, and the masses of an electron me and a baryon rriB (see, for ex- 
ample, [^). In principle, if did not not know rrie and rriB, v/e could derive 
them from accurate astronomical measurements! It seems to the author that 
in the case of relic gravitons we are dealing with the problem of a similar 
simplicity and deepness. 

Let us continue our review of Fig. 1. A spectrum, part of which is 
shown as a horizontal line at the level /i(z/) = 10~^, was derived from the 
f3 = —2 model and placed at the highest level allowed by then existing limits 



on the microwave background anisotropy AT/T. It is known ||T3[ that the 
long- wavelength cosmological perturbations produce the large- angular-scale 
anisotropy AT/T of the order of AT/T ^ h. The upper limit of that time 
AT/T < 10~^ required h < 10~^. The horizontal position of this part of the 
spectrum (use Eq. (13) for /5 = —2) explains why this spectrum is called 
"flat" (or Harrison-Zeldovich) spectrum. All waves with present frequencies 
v > vh were ordered long ago in the "fiat" spectrum, but adiabatically 
decreased their amplitudes by now. In particular, ii (3 = —2 and v > Vm 



one has hiy) ~ v^^ |14|. The full present-day spectrum is a continuation of 
the horizontal line to higher frequencies, as shown in Fig. 1. As was already 
emphasized above, we ignore oscillations in the power spectrum P{n) of the h- 
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field itself. As for the power spectrum of the energy density eg, it is expected 
to be smooth, because the sin^(n?7 + x) oscillations in /i^ combine with the 
cos^ (nr] + x) oscillations in {h'y to produce a smooth function for the sum. 
The f3 = —2 model and the AT/T limits of that time did not allow Qg in 
the frequency bands Au ^ z/ to be larger than 10~^^ for all ground-based 
and space techniques. (The old graph for the space interferometer sensitivity 
does still appear promising for detecting such a signal but it is now known 
to be overly optimistic [jl|.) 

The situation has considerably changed in the recent years. 

First, the large-angular-scale anisotropy has been actually detected [Q], 
so we are now dealing with the detected signal AT/T ^ 10~^, and not with 
the upper limit AT/T < 10~^. A very important question is whether we can 
attribute, say, a half of the detected signal to gravitational waves, assuming 
of course that the AT/T is caused by cosmological perturbations of quantum- 
mechanical origin, as we argued in the very beginning of this paper. Without 
being able presently to answer this question observationally, we should rely 
on the theory. The theory definitely says yes. A possible contribution of 
quantum mechanically generated density perturbations can be of the same 
order of magnitude as (in fact, somewhat smaller than) the gravity- wave 
contribution, but cannot be much higher [||. Specifically for models with the 
scale factors (11) governed by a scalar field, the characteristic amplitude of 
the long- wavelengths metric perturbations /i(z/) associated with the density 
perturbations and responsible for AT/T is described by exactly the same 
formula as formula (13) for gravitational waves. This is not surprising since 
the basic dynamical equation for the scalar field density perturbations 
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/i„ + [n ^^]/in = 

(where 7 = 1 + {a /a')' and the scalar field potential is arbitrary) is not 
only similar, but is exactly the same as Eq. (4) for gravitational waves, 
when the scale factor a{r]) at the initial stage of expansion is taken as a 
power-law function (11) (in which case, 7 = const, and 7 drops out of the 
equation). The numerical coefficients in formula (13) are somewhat in favor 
of gravitational waves, but it is not numerical coefficients of order 1 that we 
are now discussing. What is important for us is that the observed AT/T 
can now be taken as an experimental point for relic gravitational waves. If 
we have still used the j3 = —2 model, the entire inflationary graph of Fig. 
1 should have been shifted down by 1 order of magnitude in terms of /i(z/), 
predicting a hopelessly small Qg ~ 10~^^ for frequencies of our interest, where 
direct measurements of the gravity-wave background are possible. 

Second, the processing of the COBE data has allowed us to obtain some 
information |l^, about the power-law spectral index of primordial pertur- 
bations. Usually, the COBE data are processed under the assumption that 
the anisotropics are caused by density perturbations. Fortunately, it does not 
matter what one thinks about the perturbations while processing the data: 
the effects of gravitational waves and density perturbations on the large- 
angular-scale anisotropics are about the same when amplitudes and spectral 
indexes (/3-I-2 in our notation) of the corresponding metric perturbations are 
the same, see Eq. (13). 

Density perturbations with the Harrison- Zeldovich spectrum are often 
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desribed by the spectral index n = 1 (this n should not be confused with the 
dimensionless wave number n used in this paper). The relationship between 
this spectral index n and the parameter (3 is 

n = 2p + 5. (16) 

Information about the primordial spectral index extracted from the large- 
angular-scale anisotropics is information about gravitational waves (when 
they dominate or provide, at least, a half of the signal) even if one thinks that 
the processed anisotropics are entirely produced by density perturbations. 

The authors of concluded: "...we find a power-law spectral index of 
n = 1.2±0.3 ..." , "The power spectrum of the COBE DMR data is consistent 



with [a Peebles-Harrison-Zeldovich n = 1 universe]". The authors of [|1| 
who processed the same set of data but in a different manner, concluded: 
"The spectral parameter of the power spectrum of primordial perturbations 
n = 1.84 ± 0.29 [is] estimated", "The power spectrum estimation results are 
inconsistent with the Harrison- Zeldovich n = 1 model with the confidence 
99 %". 

It is difficult for us to judge whether the n = 1 model is ruled out at the 
confidence level 99 %, according to [|I6|, or at the confidence level 60 % or so. 



according to [T^. However, in these two results, we see the indication that 



the spectral index n is indeed larger than 1. As a compromise, we will first 
take the value n = 1.4 and derive its consequences for gravitational waves. 

The spectral index n = 1.4 translates into f3 = —1.8, see Eq. (16). The 
u < Ufj part of the spectrum is not any longer "flat" but gives more power 
to higher frequencies, hlu) ~ see Eq. (13). The value h^vn) = 10~^ 
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is fixed by observations. The number (3 = —1.8 should also be used in Eq. 
(14) and Eq. (15). This gives hiym) = 10"^-^, a slightly smaller amplitude 
than h{vrn) = 10~* for the j3 = —2 model. But for the higher frequencies 
the results are very encouraging. At the LISA-tested frequency u = 10^^ Hz 
we get h = 10-19 and Qg = IQ-^ At the LIGO/VIRGO-tested frequency 
u = 10^ Hz we get h = 10"^^ and fig = 10"®. To plot the full spectrum, 
one can essentially shift down the entire graph of the f3 = —2 model and 
slightly rotate the graph around the fixed observational point h = 10~^ at 
V = vh- This produces a solid line shown in Fig. 1. This line is no longer 
a result of a simple model-building, it is now supported, at least partially, 
by observations. The fact that the graph consists of straight line pieces 
meeting at corners is accounted for by the nature of our approximation: 
strictly power-law scale factors joined at the transition points between the 
initial, radiation-dominated, and matter-dominated eras. At a more accurate 
graph, the corners will be rounded and lines will be slightly bent. 

We have to admit that it is not so easy to give a "microphysical" ex- 
planation to the derived parameter /? = —1.8. This value of (3 corresponds 
to the equation of state p = —1.2 e, see Eq. (12). With the equation of 
state of this kind, energy density of the matter, as well as curvature of the 
space-time, increase in course of expansion. However, the value (3 = —1.8 is 
marginally consistent with the assumption that the Planck densities are not 
encountered in course of the evolution . The scalar fields, often considered 
in the context of the inflationary hypothesis, are only capable of providing 
an equation of state with f3 < —2. Possibly, a solution to the "microphys- 
ical" side of the problem can be found along the lines of the "superstring 
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motivated" cosmologies [|T^. 

The direct detection of the useful noise (a stochastic gravity- wave signal) 
by noisy detectors can be achieved with the help of a standard technique of 



cross-correlating the outputs of two or more detectors [|19], |21|. The LISA 
is not planned to have two independent detectors ^ . However, the predicted 
signal is so high, ilg = 10'^ at u = 10"^ Hz and fig = 5 x lO"*^ at z/ = 10"^ 
Hz, that one can probably recognize the signal by comparing the observa- 
tional data with the calculated sensitivity of the instrument. Fortunately, at 
frequencies around 10~^ Hz and higher, the contaminating gravity-wave noise 
from compact binaries is expected to be below the projected LISA sensitivity. 
It is important to note that the conservative value n = 1.2 (/3 = —1.9) of the 
spectral index does still lead to a quite well measurable signal: h = 10~^°'^, 
Qg = 10-^^ at 1/ = 10-3 Hz and h = 10~^^, Qg = 10"^° at z/ = 10^ Hz. 

The first evaluation |jl9| of a possibility to detect relic gravitons by a 



cross-correlating technique was based on the assumptions that the flux den- 
sity behaves as ~ i>~^ (that is, /i(z/) ~ z/"^, like in the f3 = —2 model 
considered above), that the expected amplitude is at the level (in current 
notation) of Qg = 10~^, and that the electromagnetic detectors operating in 
the high-frequency band u = 10^ Hz are being used. At that time, these 
assumptions were of kind of a stretched imagination. Presently, this pos- 
sibility may turn out to be more realistic in view of the fact that fig(z/) is 
growing toward the higher frequencies, see Eq. (15) and Fig. 1. However, the 
prospect for the high-frequency techniques, such as bars and electromagnetic 
detectors, may be not as good as it looks on the graph. The high-frequency 
parts of spectra generated in simple models of the very early Universe have 
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tendency to deviate down from the straight hne corresponding to a strict 
power-law dependence (11). Generally speaking, the farther we extrapolate 
the spectrum from the experimental point at i/ — uh, the less confident we 
are. In this respect, the LISA has an additional advantage of operating at 

relatively low frequencies. 

So, what is my conclusion? Urgently fly LISA! 
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